大 型 往复 压缩 机 管道 系统 气流 脉动 的 计算 
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摘 要 : 针对 某 大 型 往复 活塞 压缩 机 管 系 结构 ， 进 行 了 气流 脉动 的 计算 与 分 析 。 以 平面 波动 理论 为 基础 ， 采 用 网 络 传递 矩阵 方 
法 ， 建 立 了 该 复杂 管 系 气流 脉动 的 网 络 分 析 模 型 ， 并 开发 了 相应 的 分 析 软 件 程序 代码 。 计 算 了 该 管道 系统 气 柱 的 前 十 阶 固 有 频 
率 ， 并 判断 管 系 是 否 共振 ， 同 时 ， 计 算 了 管道 系统 中 各 节点 的 脉动 压力 ， 计 算 结果 满足 API 618 标准 的 要 求 。 复 杂 管 系 气流 及 
动 的 研究 ， 将 为 管 系 的 设计 、 和 气流 脉动 与 管道 振动 的 拟 制 与 消除 提供 理论 依据 。 
关键 词 : 往复 压缩 机 ， 平 面 波动 理论 ;网 络 传递 矩阵 ， 气 柱 固 有 频率 ， 脉 动 压 力 
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Abstract: The calculation and analysis of gas pulsation are carried out for a large reciprocating compressor piping 
system. Based on the theory of plane wave, a network analytical model of gas pulsation for the complex piping 
system of one large reciprocating compressor is established using the network transfer matrix method, and the 
corresponding analytical program code is developed. The first ten order of gas natural frequencies are calculated, and 
the estimate is made about the resonant situation of piping system. At the same time, the pulsating pressure at each 
calculating node of the pipeline system is calculated. The calculation results demonstrate that the requirements in API 
618 for the gas pulsating in the piping system are met. This research would provide the theoretical basis not only for 
the design of the piping system, but also for the control of the gas pulsation as well as the pipeline vibration. 
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0 引言 
大 型 往复 活塞 式 压缩 机 作为 压缩 和 输送 介质 生产 安全 造成 很 大 的 威胁 。 因 此 ， 研 究 气流 脉动 的 
的 动力 源 ， 广 泛 应 用 于 石油 、 化 工 、 制 冷 等 领域 。 产生 机 理 ， 计 算 气 流 脉 动 声场 分 布 ， 无 论 是 对 主动 


由 于 运输 尺寸 、 重 量 和 占 地 面积 的 限制 ， 对 大 型 往 宝 制 和 抑制 压缩 机 组 的 气流 脉动 ， 以 及 减 小 或 消除 
复 压缩 机 的 结构 紧凑 性 、 可 靠 性 和 运行 稳定 性 的 要 蔬 道 振动 ， 还 是 提高 运行 的 经 济 性 ， 均 具有 重要 的 
求 越 来 越 高 。 但 是 ， 由 于 压缩 机 间 短 性 的 吸 、 排 气 ， 工程 实用 价值 和 现实 意义 。 

会 使 管道 内 产生 较 强烈 的 气流 脉动 ， 气 流 脉 动 不 仅 近年 来 国内 外 对 于 气流 脉动 的 研究 ， 主 要 采用 
会 影响 压缩 机 的 工作 性 能 ， 还 会 引起 管道 振动 ， 对 转移 矩阵 法 上 四、 有 限 元 法 上 、 声 电 模拟 法 9、CFD 
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法 外 和 一 维 时 域 流动 求解 法 外 以 及 商业 软件 Bentley 略 去 平均 气流 速度 ， 可 以 得 到 平面 波动 方程 ” 


PULS。PULS 是 国外 专门 用 于 管 路 气 柱 分 析 的 软 Op, _ 1 др, E 
件 ， 目 前 已 被 国外 市 场 鞋 断 ， 不 仅 购 买 价格 昂贵 ， or ax? 
而 且 软 件 的 更 新 也 花费 巨大 ， 因 此 ， 开 发 自主 气流 其 解 为 
脉动 分 析 软 件 具 有 极其 重要 的 意义 。 转 移 矩 阵 法 ， конна ведкнаа 5 
常 称 为 网 络 传递 矩阵 方法 ， 是 目前 计算 压缩 机 管 
道 气流 脉动 最 为 实用 的 一 种 适 于 工程 设计 与 分 析 的 在 平均 气流 速度 m =0 的 情况 下 ， 根 据 运动 方 
计算 方法 。 西 安 交通 大 学 对 此 方面 的 研究 有 具有 重大 ”程式 (3〉 和 波动 方程 式 (4) 
贡献 趾 ， 特 别 是 在 气 柱 固有 频率 、 气 流 脉 动 、 压 力 e ЗЕ. (3) 
脉动 的 计算 方面 ， 成 功 解决 了 一 些 实际 问题 。 李 志 Oo 内 tr 
博 中 通过 转移 矩阵 法 ， 获 得 任意 复杂 管 系 的 固有 频 Dapp та па 
率 及 该 管 系 任意 位 置 的 脉动 波形 和 脉动 幅 值 ， 通 过 Of Ox 
比较 计算 结果 与 实测 值 , 发现 该 方法 有 较 高 的 精度 。 可 得 ， 
徐 斌 等 人 对 级 间 管 路 气流 脉动 做 了 深入 细致 的 研究 ee 
TO 中， 并 针对 传递 甜 阵 法 模型 计算 精度 不 足 的 问题 对 или? ве | > 
二 模型 进行 了 修正 加， 并 进行 实际 测量 ， 十 分 具有 指 
导 意 义 。 攀 民 外 基于 传递 箱 阵 法 计算 出 气 柱 固 有 频 式 中 ，p; ， 脉 动 压力 ，Pa ; wu， 脉动 速度 ， 


率 ， 研 究 了 管道 气流 脉动 的 控制 方法 ， 有 效 降低 了 m/s; P， 气 体 的 密度 , kg/m ; a ,气体 声速 , m/s; 
O 压缩 机 组 出 口 管道 和 相应 附件 的 振动 ， 保 证 机 组 整 A 、B"， 复 数 常数 ， 由 管道 的 端点 条 件 确定 。 
二 一” ” 体 安全 稳定 运行 ， 意 义 重大 。 复杂 的 管道 系统 可 以 看 成 是 由 等 截面 管 、 容 
下 本 文 基于 平面 波动 理论 ， 采 用 网 络 传递 矩阵 方 器 、 汇 流 点 、 异 径 管 等 基本 管道 元 件 组 成 ， 根 据 网 
法 ， 编 制 了 气流 脉动 分 析 程 序 ， 对 某 企 业 的 大 型 往 络 传输 矩阵， 可 以 建立 各 元 件 之 间 节 点 处 的 脉动 压 


复式 压缩 机 管道 系统 的 气流 脉动 进行 了 分 析 计算 ， 力 和 脉动 速度 之 间 的 关系 。 假 设 压缩 机 管道 系统 起 
得 到 管 系 气 柱 的 各 阶 固有 频率 ， 并 判断 管 系 是 否 共 。 始 端 和 末端 的 脉动 压力 和 脉动 速度 分 别 为 p,、u、 


ZN 


振 ， 以 及 网 络 元 件 各 节点 处 的 脉动 压力 ， 计 算 结果 Par u WAEHERE (6). 


S 


与 API618 标准 进行 了 比较 ， 可 以 判断 管 系 的 设计 -am (6) 

是 否 合理 。 为 压缩 机 组 管 系 的 设计 和 振动 抑制 ， 提 и ‘Lu 

供 了 理论 依据 和 改进 的 方向 。 式 中 ，M ,为 各 管道 元 件 的 传输 和 矩阵， 根据 边 
ee 界 条 件 即 可 得 到 管 系 气 柱 的 固有 频率 方程 ， 利 用 

1 数学 模型 Matlab 编程 即 可 算出 管道 系统 的 各 阶 气 柱 固有 频 
对 于 气流 脉动 的 研究 主要 分 为 两 部 分 : 计算 气 。”” 率 。 

柱 的 固有 频率 ， 使 其 避 开 压缩 机 激 振 频率 的 共振 频 1. 2 脉动 压力 计算 数学 模型 

带 (0.8~1.2) f 计算 气流 的 脉动 压力 ， 使 其 控制 在 为 计算 压力 脉动 量 ， 必 须 考虑 阻尼 。 

标准 范围 内 。 本 文采 用 网 络 传递 矩阵 法 计算 气 柱 的 线性 阻尼 波动 方程 为 : 

固有 频率 和 气流 的 脉动 压力 。 Фр Кб, ТОР 0 G5 

1.1 气 柱 固有 频率 计算 数学 模型 ar ад «а 
正常 情况 下 , 管道 中 气流 压力 的 脉动 值 相 对 于 式 中 ，R ， 管 道 沿 程 的 阻尼 系数 ， 其 表达 式 为 


Wi 4 ， 气 体 与 管 壁 间 的 摩擦 系数 ， DD ， 
道内 径 , m; ис» 管道 气流 的 平均 气流 速度 , m/s。 


均 压 的 平均 值 是 一 个 很 小 的 数值 〈 按 双 振幅 计 ， 一 
般 在 8 多 以 内 )， 符 合 平面 波动 理论 的 假设 ， 因 此 ， 


е 


往复 式 压缩 机 管道 内 的 气流 脉动 可 以 采用 平面 波 方程 (7) 的 解 为 : 
动 理 论 来 分 析 。 pi 2 [areen + Ве ®t) et (8) 


根据 气体 的 连续 方程 、 运 动 方程 和 波动 方程 ， 式 中 , w ,脉动 的 圆 频率 ; а = К/(2а) ;大 = w/a o 
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首先 对 压缩 机 的 


吸 、 排 气 曲线 进行 谱 量 分 析 ， 


得 到 各 阶 脉 动 质量 流 


量 ， 然 后 ， 根 据 各 元 件 的 传输 


和 矩阵 得 到 管 系 总 的 传 
力 ， 最 后 ， 将 各 阶 脉 
内 任 一 点 的 脉动 压力 。 
可 以 参考 文献 [1]。 


输 矩 阵 ， 进 而 求 出 各 阶 脉动 压 
动 压力 加 以 县 加 ， 可 得 到 管 系 
各 通用 元 件 的 网 络 传递 矩阵 


2 气流 脉动 计算 与 分 析 


2.1 网 络 计 算 模型 
пие т 
道 系统 内 的 介 


作用 。 
量 为 17.3 реми 


16 с, HAW 


Доли Cn 
质 是 混合 工艺 气 ， 平 均 分 子 


压缩 机 转速 为 994 r/min， 进 气压 
力 为 2 MPa ， 排 气压 力 为 6 MPa ， 进 气温 
度 为 93.9 С. 


度 为 


计算 和 建 模 中 ， 


统 作 如 下 假定 


(1) 忽略 系统 中 的 测 温和 测 压 装置 ; 
(2) 将 阀门 看 作 变 截面 ， 
缓冲 缸 看 作 体积 元 件 ; 
(3) 冷凝 器 简化 为 由 两 个 前 端 管 
箱 和 两 根 连接 管 五 部 分 组 成 ， 前 、 后 端 管 箱 看 作 体 
只 元 件 ， 连 接管 看 作 等 截面 管 ; 
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首先 对 该 压缩 机 组 排 气管 道 系 


过 滤器 、 洗 涤 饶 和 


ня. ква 


(4) 不 考虑 法 兰 ， 常 开 阀 门 均 看 作 通 管 ; 


(5) 不 考虑 排污 系统 以 及 放空 管 路 对 气 柱 的 


影响 。 


按照 以 上 假定 ， 对 实际 管道 系统 进行 了 简化 ， 
1 所 示 ， 即 为 分 析 管 系 气流 脉动 的 网 


简化 结果 如 图 
络 阻抗 模型 。 
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个 单元 ，46 个 节点 ， 


16 分 别 为 六 个 气缸 的 排 气 端 ， 节 点 40 为 声学 开口 


边界 条 件 〈 工 艺 气 出 
条 件 (氮气 置换 口 ， 
门 关闭 )。 


22 气 柱 固有 频率 的 i 


在 气流 脉动 分 析 时 


图 1 管 路 系统 网 络 模型 


Fig.1 Network model of pipeline system 


该 管 路 系统 被 划分 为 45 
其 中 节点 1、3、5、12、14、 


口 )， 节 点 46 为 声学 闭口 边界 
一 般 情况 下 该 分 支管 路 上 的 阀 


计算 


对 于 声学 上 的 开口 有 脉动 压力 p=0 ,声学 上 的 


闭口 有 脉动 速度 u =0，, 把 端点 不 为 零 的 p 或 u 设 为 
1， 与 气缸 相连 的 管 路 一 端 


看 成 是 闭 端 。 将 0.8~1.2 


倍 气 柱 固有 频率 之 间 


的 前 十 阶 固 有 频率 f 
1 所 示 。 


根据 传输 矩阵 法 计算 出 排 气 管道 


的 区 域 定义 为 共振 区 。 
首 系统 内 气 柱 
及 共振 区 域 (0.8 ~1.2) 7, ПА 


表 1 气 柱 固有 频率 及 共振 区 


Table 1 Gas natural frequency and resonance region 


n f/ĦHz (0.8~1.2)f/ Hz 
1 0.6958 0.56~0.83 
2 2.3586 1.89~2.83 
3 3.1020 2.48~3.72 
4 5.1612 4.13~6.19 
5 7.9449 6.36~9.53 
6 9.9331 7.95~11.92 
7 12.6257 10.10-15.15 
8 15.3800 12.30-18.46 
9 17.8260 14.26-21.39 
10 19.1640 15.33-23.00 


压缩 机 的 激发 频率 уа (9) 计算 ， 匀 度 应 限制 在 7% 以 下 或 小 于 下 面 公式 的 计算 值 。 


fs на. (9) [8] =0.03R x100% (11) 
式 中 ，N ， 压 缩 机 转速 ，r/min; m ， 单 作用 时 式 中 ，R ， 级 压力 比 。 
等 于 1， 双 作用 时 等 于 (b) 缓冲 装置 管 侧 接管 处 或 接管 以 外 的 最 大 
НТ 的 压缩 机 的 主 激发 频率 为 33.13 Hz， 许 用 脉动 范围 
通过 计算 结果 可 以 看 出 ， 压 缩 机 的 激发 频率 有 效 的 缓冲 装置 应 将 缓冲 装置 管 侧 的 压力 不 均匀 度 
避 开 了 气 柱 前 十 阶 共振 区 ， 说 明 不 会 引起 压缩 机 与 限制 在 式 〈12) 范围 内 ， 
管道 内 的 气 柱 共振 。 ET | е 
2.3 脉动 压力 的 计算 д|- Va/350 (12) 
2.3.1 压力 脉动 的 控制 标准 
压力 脉动 的 程度 用 压力 不 均匀 度 56 来 评价 , 即 ， 式 中 ，a ， 气 体 的 音速 ，m/s; P ， 管 路 平均 
ро Виа _ 2AP do 绝对 压力 ，MPa;， D> ETW, mm; f, ВК 
R R 频率 ，Hz。 
式 中 ，P,, ， 气 体 瞬时 压力 最 大 值 ，MPa ; Ра» 2.3.2 脉动 速度 谐 量 分 析 
气体 瞬时 压力 最 小 值 ，MPa ; p,， 气 体 平 均 压 力 节点 1、3、5、12、14、16 为 压缩 机 排 气 端 边 
值 ,，MPa; Ap ， 压 力 脉动 最 大 幅 值 ，MPa 。 且 界 条 件 ， 脉 动 速 度 为 已 知 ， 首 先 对 脉动 速度 和 脉动 
O A, ВАР +Р.)/2, AP=(P, -Po,)/2， 所 质量 流量 进行 傅 里 时 谐 量 分 析 ， 节 点 1 处 的 脉动 速 
© идр-р6/2. 度 和 脉动 质量 流量 前 10 阶 谐 量 如 表 2 所 示 。 其 余 排 
= 根据 2007 年 发 布 的 第 五 版 API618， 系 统 中 各 气 端点 的 各 阶 脉 动 速度 的 谐 量 系数 不 同 ， 但 是 幅 值 
T 处 的 压力 脉动 控制 标准 概述 如 下 09。 都 与 节点 1 处 相同 ， 在 此 不 再 一 一 列 出 。 由 表 中 可 
даа Са) 压缩 机 气缸 法 兰 最 大 许 用 压力 脉动 以 看 出 ， 对 双 作 用 气缸 来 说 ， 最 大 幅 值 的 谐 量 为 第 


= 压缩 机 气缸 法 兰 处 未 经 滤波 的 压力 脉动 不 均 2 阶 ， 其 频率 就 是 激发 主 频 率 。 


表 2 节点 1 脉动 速度 和 脉动 质量 流量 前 10 阶 谐 量 


Table 2 The first 10 order of harmonics for the pulsating velocity and mass flow rate at node 1 


c 脉动 速度 分 量 脉动 质量 流量 分 量 
п 
а jms- b, /m s?! с, /m s! Go Kes Š, kg-s! č, рев“ 

1 0.2534 -0.2132 0.3312 0.1328 -0.1117 0.1736 
2 0.0284 -2.4228 2.4230 0.0149 -1.2697 1.2697 
3 -0.1230 -0.1997 0.2346 -0.0645 -0.1047 0.1229 
4 -1.4266 0.1531 1.4348 -0.7476 0.0803 0.7519 
5 -0.0867 0.0530 0.1016 -0.0455 0.0278 0.0533 
6 0.4410 0.4487 0.6291 0.2311 0.2351 0.3297 
7 0.0268 -0.0075 0.0278 0.0140 -0.0039 0.0146 
8 -0.1369 -0.5336 0.5509 -0.0717 -0.2796 0.2887 
9 -0.0391 -0.0172 0.0427 -0.0205 -0.0090 0.0224 


10 -0.3073 0.3382 0.4570 -0.1610 0.1772 0.2395 
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图 2 为 一 个 周期 内 脉动 速度 的 曲线 图 。 从 周期 x 104 
内 脉动 速度 曲线 图 可 以 看 出 该 图 形 符合 双 作用 气 饶 
的 排 气 速度 规律 。 同 时 ， 也 可 以 看 出 在 压缩 机 的 i 
侧 和 盖 侧 的 脉动 速度 幅 值 是 不 对 称 的 ， 这 是 由 于 轴 а 0 
侧 活塞 杆 占据 了 一 定 的 气缸 容积 ， 结 构 不 对 称 导致 pe 
的 。 
-10 
0 1 2 3 4 5 6 
0/гад 
j 图 3 压缩 机 排 气 端点 1 的 压力 波形 图 
3 4 Fig.3 Pressure waveform of node 1 at exhaust pore 
Na 
> 2 表 4 EHO A 7 的 压力 幅 值 
0 Table 4 Pressure amplitude of node 7 at Ше outlet of Ше buffer 
0 1 2 3 4 5 6 n Re/Pa Im/Pa Pm/Pa d/rad 
Wrad 1 176.18 587.8 613.64 1.2796 
2 脉动 速度 曲线 2 1504 9369.5 9489.4 1.4116 
Fig.2 Pulsating velocity сигуе 3 184.45 1567.7 1578.5 1.4537 
233 压力 不 均匀 度 的 计算 4 1606.8 14402 14491 1.4597 
根据 各 元 件 的 单元 传输 矩阵 就 可 以 得 到 各 点 3 
各 阶 脉动 压力 谐 量 ， 由 于 节点 太 多 ， 在 此 只 列 出 压 сл оо а a 
GEHE НАСОКИ AEJ, WE 3 E ли со ареста 
8 55221 -46916 72460 -0.7043 


FIK 4 所 示 。 


将 各 阶 脉动 压力 幅 值 按 它们 的 相位 关系 登 加 
起 来 ， 得 出 端点 7 总 的 压力 波形 图 ， 见 图 4。 其 中 


表 3 压缩 机 排 气 端点 1 的 压力 幅 值 


Table 3 Pressure amplitude of node 1 at exhaust pore 


脉动 压力 最 大 值 为 1564.6 Pa， 最 小 值 为 -1200 Ра, 


пъ O ”一 其 中 最 大 最 小 脉动 压力 分 别 出 现 在 转角 071 rad 
1 1.92 125.45 125.47 1.5555 和 5TA tad a 管内 平均 压力 为 6 MPa, 其 压力 不 
2 1.67 700.35 700.54 1.5470 
3 0.97 25.65 25.67 1.5328 SIEA 0.046 %。 
4 2.66 95.15 95.19 1.5428 
5 0.94 10.88 10.92 1.4844 1000 
6 -36.94 -255.35 258.01. -1.7145 B00 
7 -2.15 -6.34 6.70 -1.8976 ` o. 
8 392.05 -300.44 493.92 -0.6539 
EP, Re RUR: Im RER, Ра EIIE: Дон 
d， 相 位 角 。 一 一 ста 


3 
Orad 
4 ”缓冲 钠 出 口 点 7 的 压力 波形 图 


Fig.4 Pressure waveform of node 7 at the outlet of the buffer 


将 各 阶 脉动 压力 幅 值 按 它们 的 相位 关系 车 加 
起 来 ， 得 出 端点 1 总 的 压力 波形 图 ， 见 图 3。 其中， 
脉动 压力 最 大 值 为 81700 Pa， 最 小 值 为 -102910 Ра, 
最 大 最 小 脉动 压力 分 别 出 现 在 转角 2.26 rad 和 5.82 依次 类 推 ， 可 以 求 出 管道 中 各 点 脉动 压力 的 最 
rad 处 。 管 内 平均 压力 为 6 MPa,， 其 压力 不 均匀 度 为 大 值 和 最 小 值 ， 进 而 求 出 压力 不 均匀 度 ， 并 与 
3.08 %。 API618 标准 进行 比较 ， 如 表 5 所 示 。 


paj 


表 5 各 节点 的 压力 不 均匀 度 


Table 5 Pressure un-uniformity of each node 


节点 入 点 位 置 已 ， /Pa P” „Ра 已 /MPa 6 /% д,/9 

1 气 饶 排 气 端 口 81700 -102910 6 3.08 7 

3 气 饶 排 气 端 口 82030 -100800 6 3.05 7 

5 气缸 排 气 端 85100 99990 6 3.08 7 

7 Zh аи 1564.6 -1200 6 0.046 0.51 
12 气 饶 排 气 端 97850 -88050 6 3.10 7 

14 气 饶 排 气 端 口 81700 -98070 6 3.00 7 

16 气 饶 排 气 端 口 82820 -101250 6 3.07 7 

18 Сари fan 1234 -982 6 0.037 0.51 
表 中 ，P”",,、， 脉 动 压 力 幅 值 最 大 值 ， P”,,，， 脉 动 压 力 幅 值 最 小 值 ， 号 ， 脉 动 压力 平均 值 ，6 ， 压 力 不 均 匀 度 ;6,， 


API 618 允许 值 。 


从 表 中 可 以 看 出 ， 经 过 缓冲 鳞 后 压力 不 均匀 度 号 :0214120042) ,特此 感谢 ! 
远 小 于 API618 允许 值 ， 因 此 ， 后 面 的 管道 中 压力 
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